Quantum States of Matter of Simple Bosonic Systems: 
BEC's, Superfluids and Quantum Solids 
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The phase diagram of a single component Bose system in a lattice at zero temperature is obtained. 
We calculate the variational energies for the Mott insulating and superfluid phases. Below a certain 
critical density the Mott insulating phase is stable over the superfluid phase for low enough tunelling 
amplitude regardless of whether the number of bosons is or is not incommensurate with the lattice. 
The transition is discontinuous as the superfluid order parameter jumps from a finite value to zero 
at the Mott transition. 
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The purpose of the present work is to study the phase 
diagram of a single component Bose gas at zero temper- 
ature. Specifically, we have in mind a system of bosons 
in an optical lattice, where the wells on each site are 
made deeper thus causing a transition from the super- 
fluid phase into a phase where bosons are localized on 
the lattice sites, i.e a quantum solid or Mott insulating 
phasei*^. Recently such a transition has been observed 
exp er iment allyi. 

The original theoretical work by M.P.A. Fisher et. 
al^ on the superfluid-Mott transition in bosonic sys- 
tems predated the experimental realization of artifi- 
cially engineered Bose-Einstein condensate o'^'^'^ i(BEC's), 
which are superfluid systems as they posses a fi- 
nite sound speedSiii*i4 and suppressed long-wavelength 
scattering^. In that original work the transition was 
asserted to be continuous and possible only for a num- 
ber of bosons commensurate with the lattice, i.e. in- 
tegral number of bosons per site. While widely be- 
lieved to be true, both assertions are quite puzzling and 
may contradict experiments^. In fact, recent theoretical^ 
and numerical workiSi has pointed to a transition that 
is incommensurate and "nonstandard" even for a pure 
systemifl. In the present note we extend the work in the 
pure system within the Bogolyubov approximation for 
the superfluidii*i^ 

The superfluid-Mott transition is between a superfluid 
and a quantum solid. In systems without disorder, fluid- 
solid transitions are usually discontinuous as the order 
parameters are too differenliii. In the superfluid system 
that is certainly the case. The Mott insulating phase 
is characterized by a well defined number of bosons per 
site, one could take the density as its order parameter. 
The superfluid phase has off-diagonal order with break- 
ing of gauge invariancei^, i.e. U{1), characterized by a 
coherent ground state of Bogolyubov pairsii. The tran- 
sition will be characterized by a discontinuous jump in 
the superfluid response^^. 

We calculate the variational energies of Mott insulat- 
ing and superfluid phases with both commensurate and 
incommensurate number of bosons. We find that below 
some critical value Vc of the ratio t/U the Mott insulating 



phase is energetically favorable over the superfluid phase 
as long as the number of bosons per site is not too high, 
i.e. less than a critical value ri^. Below and the 
superfluid phase is never stable regardless of commensu- 
rability. 
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FIG. 1; Schematic Phase Diagram of a single component 
Bose system at zero temperature, n is the average number 
of bosons per site. 



When the number of bosons per site is higher than the 
critical number Uc, the quantum solid phase is never sta- 
ble within our approximations and there is only a super- 
fluid phase even for arbitrarily small t/U. On the other 
hand, we believe our approximations, which are apt at 
low number of bosons per site, fail at high densities and 
there will be "crossover" to physics that is more like that 
of Josephson coupled arrays:^' This is a matter for fur- 
ther study. The phase diagram is shown in figure 1. 

A very apt prototype Hamiltonian to study bosons in 
a lattice (which, of course can be an optical lattice) is 
the Bose Hubbard Hamiltonian: 
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H = ^ {alaj + a]ai) + C/^ a|aj(a|a^ - 1) 



(1) 



where t is the nearest-neighbor tunneUng or hopping am- 
phtude, U is the on-site repulsion, and e is the well depth 
for each lattice site. We will suppose our lattice to have 
N sites and our system to have M bosons with M/N no 
necessarily an integer. 

When U — 0, the ground state of such a Hamiltonian 
is a BEC with all atoms in the zero momentum state. 
Such a system is not superfluid as its excitation spectrum 
has quadratic dispersioni^. An arbitrary small amount 
of repulsion causes the system to order and become su- 
perfluid. The ground state, first found by Bogolyubo\iii, 
is 
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where 



with elk = Ek + ^^U, where Mq is the number of parti- 
cles in the condensate, and = — 2t(cos k^a + cos kya + 
cos kza)+6t is the dispersion in the lattice which goes like 
ta^k^ = h^k"^ /2m in the long wavelength limit. In order 
to perform our calculations we will introduce a cut-off 
defined by the point at which the kinetic energy becomes 

comparable to the on-site repulsion, i.e. fee = yf^^" ■ 
In the Bogolyubov state, "conservation" of bosons im- 
poses the sum rule 
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which determines the number of particles in the conden- 
sate. In the Bogolyubov state we have a superfluid order 
parameter which corresponds to boson pair correlations 



with _Ek = \/(^y fik + 4:^11 the quasiparticle excitation 
spectrum of the superfluid. Since at long-wavelengths 
it is a sound spectrum, the system superflows by the 
Landau argumenfei^. Even though a Bose condensate is 
necessary for the Bose system to be a quantum liquid, 
the Bogolyubov pairs arising from the repulsion between 
atoms give the system the necessary rigidity to be su- 
perfluid. Therefore the superfluid order parameter is the 
amplitude of Bogolyubov pairs as whenever this ampli- 
tude vanishes, the system dissipates. 
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FIG. 2: (a) Dependence of the order parameter on the fraction 
of bosons in the condensate, (b) U/t vs. Mo/M 

In BEC's in an optical lattice the Mott transition is in- 
duced by making the wells on the lattice sites deeper (e 
is the lattice depth) thus suppressing the tunelling am- 
plitude t between near neighbor sites, i.e. making the 
bosons "heavy" . Thus the relevant ratio of on-site repul- 
sion U to hopping t increases. As this ratio increases, 
the superfluid order parameter increases as shown in the 
figure 2. The increase of the order parameter means that 
the superfluid phase is becoming more stable over the 
Bose condensed gas. On the other hand, we also plot 
the number of particles in the condensate. This number 
is becoming small as U/t increases. As the condensate 
gets depleted the system is becoming "less fluid" and will 
solidify into a Mott insulating phase with a well deflned 
number of bosons per site. Since the Bogolyubov order 
parameter will be zero in the quantum solid phase, there 
will be a discontinuous jump in the superfluid response 
at the Mott transition. 

The Mott-insulating ground state wavefunction for 
pN < M < {p + l)N, with p, an integer is 
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where and aj are boson creation operators at lattice 
(5) sites ij and I respectively. For an incommensurate sys- 
tem, pN < M < {p + 1)N, the ground state will not be 
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strictly periodic as there are M — pN that do not fill all 
lattice sites. No two of these extra bosons would occupy 
the same site when the system is in its ground state, but 
otherwise they will be randomly distributed. In real life, 
besides on-site interactions, there are interactions acting 
at long distances. They can be attractive or repulsive. 
These interactions are weaker than the on site interac- 
tion U, hence not important for a commensurate system. 
But if there is incommensuration, these interactions will 
probably play a role in determining the positions of the 
extra bosons. 

Let us take a look at the expectation value of the su- 
perfluid order parameter in the Mott phase: 

<v> 

where 

{0\l[ara.a,l[ialr\0)= I] (Ol^^.^f !«) («) 

I I 

In the expression above, qi indicates the number of 
bosons in site I and Fij is a function of bosonic oper- 
ators on sites i and j. 



insulating phase is to have a discontinuous change in the 
superfluid order parameter. The Mott transition is dis- 
continuous, i.e. first order, with the superfluid order pa- 
rameter, and thus the superfluid response, experiencing 
a sudden jump to 0. 

The transition into the Mott insulating phase will oc- 
cur when the Mott wavefunction © has less variational 
energy than the Bogolyubov superfluid wavefunction ^ 
With the Hamiltonian given by the energy of the 
Mott state is found to be 



i,S 

-{e + U)J2{^l\n^\i^l)+uJ2{Hn^\i^^) (12) 

i i 

The first term corresponds to the kinetic energy. It de- 
stroys a boson in site i and creates it in its nearest neigh- 
bor site i + 6. Thus the lack of overlap of the resulting 
wavefunction with the ground state gives 0. We get 

{^i\n\i^i) =-{e + U)M + U{M - pN){p + 1)2 
+ U[N - (M - pN)]p^ 



+ Y[m,ar\airai\o) 

(9) 

Continuing the process we find 

- n qmaTafa,a,{alri<^]r'\0) (10) 

Since there is a non matching number of operators acting 
on the ground state this yields 0. The Bogolyubov order 
parameter is in the Mott phase 

(flka-k) = (11) 

Actually this is obvious because the Mott state has a well 
defined number of bosons per site and number phase un- 
certainty makes the off-diagonal order zero. We thus see 
that, as anticipated, the only way to go into the Mott 



(V'llHlV'i) = -sM + 2pUM -{p+ l)pNU (13) 

We see that 2pUM > {p+l)pNU. After plenty of tedious 
algebra, with our cut-off method the variational energy 
the superfluid ground state can be estimated to be 

- 6tM (14) 

For U/t small, we see that the interaction terms in the su- 
perfluid are almost irrelevant when compared with those 
of the Mott phase so the superfluid state is the stable one. 
We also notice that for not too large an M the interaction 
terms in the superfluid grow a lot faster than the interac- 
tion terms in the Mott insulating phase. Hence for some 
critical U/ 1 the Mott insulating phase will become stable 
causing the system to undergo a quantum phase transi- 
tion. Whether M is commensurate or incommensurate 
with the lattice is irrelevant to the energetics: The state 
with all bosons localized is energetically favorable to the 
state with all bosons superfluid. Whether the instabil- 
ity of the superfluid is toward a phase with all bosons 
localized, or into a state with coexistence of a commen- 
surate Mott insulator with the leftover bosons superfluid 
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is not straightforward to answer variationally. It can be 
answered experimentally as in the case of coexistence the 
superfluid response will have a discontinuous jump to a 
smaller nonzero value, while in the case of all bosons lo- 
calized the superfluid response will have a discontinuous 
jump to 0. 

The phase diagram of the simple Bose system at zero 
temperature shown in figure 1 differs considerably from 
the one proposed in the pioneering theoretical work on 
Bose system^S and thus requires some comments. The 
first and, perhaps, most important difference is that we 
do not consider a phase only model as it is not appro- 
priate to the Mott insulator at low densities. In that 
original work, the phase of the superfiuid was disordered 
by the increasing repulsion thus leading to a transition. 
In such a transition the solid would have a supefiuid order 
parameter within each well, but the phase of the order pa- 
rameter will have become randomized exactly analogous 
to what happens in Josephson junction arrays when the 
charging energy is sufficiently large. Such an insulating 
phase does not correspond to the Mott insulator we stud- 
ied here and cannot exist at small number of bosons per 
site as one needs a macroscopic number of particles to 
have a superfluid order parameter. Finally, we stress the 
the phase only model studied in the early workSii& is a 



correct model of an array of Josephson coupled superfluid 
systems and should work for bosons in an optical lattice 
when the number of bosons per site is large enough to 
make superfluid order within each well possible. Joseph- 
son coupled systems can easily be studied in an optical 
latticeiSiSi. It will be extremely interesting to see what 
the experimental phase diagram turns out to be. 

The correct physics for the Bose system phase diagram 
can be easily differentiated experimentally as the origi- 
nal work would predict that the superfluid response goes 
continuously to zero at the transition. We predict that 
the superfluid response will have a sharp jump to zero at 
the transition. We predict a Mott transition irrespective 
of commensurability as long as the number of bosons per 
site is less than some critical density Uc- At high den- 
sities there should be a "crossover" into the physics of 
coupled Josephson array a^i^^ . 
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